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We present a consistent description of the current distribution in the quantum Hall effect, based on 
two main ingredients: the location of the extended states and the distribution of the electric field. 
We show that the interaction between electrons produces a boundary line below the Fermi energy, 
which extends from source to drain. The existence of this line and that of a physical boundary 
are responsible for the formation of a band of extended states that carry the Hall current. The 
number and density of these extended states are determined by the difference between the energy 
of this equipotential boundary line and the energy of the single extended state that would exist in 
an infinite system. This is used to prove that the band of extended states is distributed through 
the bulk of the sample. We explore the distribution of the Hall currents and electric fields in by 
presenting a model that captures the main features of the charge relaxation processes. Theoretical 
predictions based on this model and on the preceding theory are used to unambiguously explain 
recent experimental findings. 
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I. INTRODUCTION 

Since the discovery of the quantum Hall effect (QHE)0 
two different approaches have been used to describe-: 
one focuses on the properties of the edge channels! 
while the othfic. concentrates on the states in the bulk 
of the systemma. Both of these approaches associate the 
quantization of the transverse conductivity o~ xy in units 
of e 2 /h, and the absence of dissipation (a xx — 0) with 
the formation of an incompressible liquid state in the 
bulk of the system, and both are based on the existence 
of delocalized states. However, the nature of the injected 
Hall current i^.quite different in these theories. In the 
first approachnd, which neglects inter-electron interac- 
tion, electrons from the current source are injected into 
delocalized states at the Fermi level which carry the Hall 
current. Since the only states at the Fermi level are lo- 
cated at the edge of the system (edge channels) , the Hall 
current has to be confined to this narrow region near the 
sample boundary. The magnitude of this edge current 
is ne 2 /h times the difference between the chemical po- 
tentials at the two edges of the system. In the latter 
approach the Hall current density is proportional (with 
the same constant of proportionality) to the electric fwld, 
whose appearance is due to the charging of the systernacl. 
In this case the Hall current is still the same but it gen- 
erally flows both along the edge and in the bulk of the 
system depending on the Hall electric field distribution. 

Experimentally, it has proved difficult to determine 
which of these descriptions is correct, and different probes 
produce apparent contradictory answers. Measurements 
of thc-equilibration rates between the current carrying 
statesEl, as well as experiments on non-local conductiv- 
itytil have natural explanations in terms of edge currents. 



However, studies of the breakdown of the QHE0, the, 
direct measurements of the Hall voltage in the sample 
and recent measurement of electric field distributionE-3 t£ 
favor the picture of bulk currents. It is the objective 
of this paper to clarify the situation and provide a de- 
scription of the current distribution in the quantum Hall 
regime which is consistent with the seemingly opposing 
experimental results. 

It was proved earlicr0 that the Hall current is absent 
in the regions near the edge where the electron density 
increases from zero to its bulk value. These are the con- 
ventionally defined edge channels. Yet one can imagine 
the situation when most or all of the current does flow 
right next to these regions, along their boundaries. More- 
over, the existence of the disorder potential modifies the 
shape and the whole notion of edge channels as they can 
no longer be thought of as simple homogeneous strips 
'parallel' to the boundaries. 

The complexity of the near-edge as well as of the bulk 
geometric structure was revealed by recent experiments 
on imaging the potentiaLdistribution in the QH bar at 
different filling factorsEjli3. It is worth noting that the 
relation between the Hall electric field and Hall current 
in such a complex inhomogeneous system can only be 
reasonable with some kind of averaging procedure. Lo- 
cal currents are the combination of non-equilibrium Hall 
current and diamagnetic currents and, in general, it is 
not possible to differentiate them without averaging over 
distances of the order of a typical localization length of 
the closed trajectories. A similar argument holds for the 
disorder and Hall voltage contributions to the electric po- 
tential. Therefore one has to be very careful inferring the 
Hall current distribution from the measured electrostatic 
potential distribution. When the Hall voltage is suffi- 
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ciently large the external electric field smears out some 
of the disorder potential fluctuations and observation of 
this electric field distribution on shorter scales becomes 
possible. However, the state of the QH system in such 
strong Hall electric field is far from equilibrium, and the 
information obtained from non-linear regime measure- 
ments, although interesting and valuable in itself, does 
not interpolate into a linear regime (experiments associ- 
ated with the study of the breakdown of the IQHE pro- 
duced current- voltage characteristics yew .distinct from 
linear at sufficiently high Hall voltageaLfio) . 

The problem of the current distribution in the QH sam- 
ple is also intrinsically connected with the question of the 
existence of delocalized states in the presence of.disarder. 
The latter was widely discussed in the literaturecjEJ but 
was not usually related to the problem of current distri- 
bution. TrugmanEU showed the connection between the 
existence of delocalized states and the classical percola- 
tion problem in an infinite-system with long-range disor- 
der potential. In his paperEil and in Ref. [20]thc possibility 
of creation of current carrying extended states by exter- 
nal electric field is discussed. The role of the physical 
boundaries of the system on the existence of delocalized 
states was also addressed in Ref. ^0|. However, it was still 
not clear where these delocalized states are located with 
respect to the physical boundaries of the system, and 
whether the derealization due to the applied external 
electric field is confined to the edges or not. 

In this paper we examine the distribution of extended 
trajectories in equilibrium and use these findings to study 
the distribution of the Hall current in the bulk. The 
conditions for the electric field to be concentrated at 
the boundaries or spread into the sample arc found and 
compared with corresponding experimental values. We 
find that, in addition to the restrictions put forth in our 
previous workEJ, the Hall current paths are distributed 
throughout the bulk of the system under the most com- 
mon experimental conditions. The interplay between the 
inter-electron interaction and the long-range disorder po- 
tential turn out to be crucial in establishing such a cur- 
rent distribution. 

The paper is organized as follows. In Sections |l| and HI 
we discuss the equilibrium state of a Quantum Hall sys- 
tem when no Hall current is present. We prove that there 
are delocalized states in incompressible regions both near 
the edges and in the bulk of the sample. Section IV 



deals with the Hall field and Hall current distributions 
away from equilibrium. We find that, similarly to the 
situation in equilibrium, extended states are distributed 
throughout the sample, their particular spatial density 
being dependent on disorder configuration and on the 
rate of relaxation processes. In Section |v| we discuss the 
validity of Stf eda formulae in the context of previous ar- 
guments. The controversy around the formula is resolved 
using the arguments of the previous sections. Finally, in 
Section VI a consistent interpretation of recent experi- 
ments is presented. Although the findings of these exper- 
iments do not appear to be in complete agreement with 



the absence of near-edge Hall current, we demonstrate 
how these results confirm the theory of the distribution 
of extended states presented in this paper when all the 
experimental conditions are fully accounted for. 



II. COMPRESSIBLE AND INCOMPRESSIBLE 
BULK 



Consider a typical system where the QHE is observed: 
a 2DEG formed at the interface of the GaAs/GaAlAs het- 
erostructure, placed in a strong magnetic field B, with 
the length of the system L y being much larger than 
its width L x . In absence of disorder, in the Landau 
gauge, with vector potential A y — B x, the electrons are 
free in the y-direction and quantized in the x-direction. 
Wave functions can be written as e lky f n (x — kl 2 ), where 
I = [Tic/ eB) 1 / 2 is the magnetic length and n enumerates 
the Landau levels (LL's). In homogeneous systems, wave 
functions and energies are simply those of the harmonic 
oscillator: e„ = huj c (n + 1/2), where uj c = eB/mc is the 
cyclotron frequency. 

The disorder potential plays a crucial role in estab- 
lishing the QHE. Two different types of disorder po- 
tential, short-range and long-range, require different ap- 
proaches to elucidate the properties of the 2DEG. The 
former is realized, e.g. in the MOSFET-devices and 
InAs heterostructures; it leads to low- mobility samples 
and is known to produce wide QH plateaus and scal- 
ing in the transitional regimes between the QH plateaus 
even at temperatures above 1 K c3. High-mobility sam- 
ples are formed at the interface of GaAs-Al^Gai-^As 
heterostructures, characterized by a long-range disorder 
potential. This potential is created by a Si-ion layer sepa- 
rated from the 2DEG by .a_undoped spacer with a typical 
thickness d ~ 100 — 500AE3. In this paper we will be con- 
cerned only with this latter type of the systems. In con- 
trast to low-mobility samples, the inter-electron interac- 
tion is very important in the case of a long-range disorder 
potential: electrons tend to screen potential fluctuations 
whose wave-length A is larger than both the magnetic 
length (see Refs. |8 29j3?]). At the same time the finite 
distance to the donor plane strongly suppresses fluctu- 
atiojasgith wavelengths shorter than the spacer thick- 

In presence of a strong magnetic field, however, screen- 
ing is limited by the incompressibility associated with the 
existence of the gap between the LL's, and filled LL's 
do not contribute to screening. Still, if the filling factor 
v = ttPtlq, no being the electron density, is not close 
to an integer, there are plenty of unoccupied states on 
the top-most occupied Landau level, and the long-range 
fluctuations of the random potential get screened in a 
significant part of the sample areac3. When the filling 
factor is close to an integer one expects to find regions 
in which the fluctuations of the random potential are not 
screened. In these regions some LL's are completely filled 
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while the rest are completely empty. When a region with 
N completely filled LL's percolates throughout the sam- 
ple the QHE is observed. The number N is usually the 
integer closest to the average filling factor v. 

Because this percolation picture of the QHE is valid 
only for filling factors close to an integer, it cannot de- 
scribe the whole range of magnetic fields where the QHE 
is observed. Evidently, some other aspects of this prob- 
lem, not accounted for in this simple picture, lead to 
the QHE at occupation numbers away from integers. 
This impoirtant issue will be discussed in an upcoming 
publicationE2l. 

In this work, we restrict ourselves to systems in which 
the region occupied by the incompressible liquid perco- 
lates, and study the current distribution in this percolat- 
ing region. The potential at every point of the sample 
is determined by the random distribution of impurities 
in the dopant layer and by the distribution of the elec- 
tron charge all over the 2DEG-plane. The potential thus 
generated has only relatively long-wavelength variations, 
and this, combined with the strong magnetic field, al- 
lows for the following simplification: electrons can be 
considered-drifting along the lines of the constant classi- 
cal energyo. The Hall current injected into the sample is 
therefore carried by the electrons moving along the tra- 
jectories stretching from source to drain. Thus, there are 
two factors that eventually determine the spatial distri- 
bution of current in the sample: the steady state charge 
and Hall electric field distributions, and the location of 
the extended (i.e. percolating from source to drain, the 
definition which we will be using throughout the paper) 
classical trajectories. To study these properties in the 
steady state, we first address the same questions for the 
equilibrium state of this 2DEG system. 



III. EXTENDED TRAJECTORIES AND 
DIAMAGNETIC CURRENT DISTRIBUTION IN 
EQUILIBRIUM 

In this section we develop a consistent description of 
the current carrying states in equilibrium. It is proven 
that extended states are present both deep in the bulk of 
the sample and close to its edges. 



A. Edge channels in presence of disorder: definition 
and absence of Hall current 



Let us consider an infinite two-dimensional random 
potential. It is clear from percolation theory that all 
equipotential lines but one are closedEJ. As a corollary 
we find that in an infinite 2DEG in a strong magnetic 
field for each LL there is a single energy Ei such that 
the state with this energy is extended in both directions. 
The only percolating equipotential line corresponds to 





FIG. 1. Finite size effects on the extended states, a. Elimi- 
nation of percolation in the longer {/-direction. Trajectory AB 
is one from the band of the states percolating in x-direction; 
line CD is part of a closed trajectory in the corresponding infi- 
nite system, b. Effects of a confining potential: lines AB and 
CD link to each other to become a single closed equipotential 
line. Other trajectories are also deformed by the confining 
potential. 



this extended state (one per each LL). If a finite-size sys- 
tem were equivalent to a piece of the same size in the in- 
finite system (with no inhomogeneities introduced by the 
boundaries), a band of states around Ei whose localiza- 
tion length is larger than the sample size, would appear 
as delocalized. However, in the finite size sample these 
states can percolate only in one direction, namely, the 
shorter oneo. The simple assumption above produces a 
drastic consequence: there are no extended states in the 
longer y-direction, and we should expect no Hall current 
in this direction (Fig. 0a). It is clear that the boundaries 
produce a non-trivial effect. 

In reality, the 2DEG is confined in the shorter re- 
direction by a self-consistent potential created by the pos- 
itive background, electron charge distribution, and gate 
voltages. Such a potential leads to an increase of the 
density of electrons from zero at the physical edge of the 
2DEG to the bulk value. Depending on the sharpness of 
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the confining, potential such growth can be either gradual 
or step-likeEjH. For illustrative purposes, let us assume 
that the charge density grows gradually (but the argu- 
ments that follow are valid in either situation). In this 
case, a narrow compressible region with a partially filled 
LL percolating in y-direction is formed near the edgeEil. 
In this region the variations of the potential due to con- 
finement and disorder are largely screened by electrons. 
As long as there are enough unfilled states on the same 
LL, electrons can screen variations of the potential. The 
density profile n(x) in the 2DEG at the edge is approxi- 
mately given by the electrostatic solutiorH 



n{x) 



I 



x + 1 



1/2 



no , x > I. 



(1) 



Here x is the distance from the physical boundary of the 
2DEG and / is the width of the depletion layer. However, 
as the occupation number of the first filled LL grows as 
required by the electrostatic solution (^), less and less 
'space' is available for variations of the local electron den- 
sity, and consequently, more and more variations of the 
potential are left unscreened. Finally, as we move to- 
wards the bulk of the system, the electrons completely 
loose the ability to screen the potential fluctuations, and 
these unscreened regions link to form a continuous strip 
of incompressible liquid along the y-direction. The out- 
ermost incompressible strip corresponds to a single filled 
LL. Moving further away from the edge, the potential 
starts falling. If more than one LL are filled in the bulk, 
the potential drops by Huj c in each incompressible strip, 
and then another compressible river (with partial occu- 
pation in the next LL), is formed roughly parallel to the 
y-directionlEl. When the last (iV-th) LL becomes com- 
pletely filled, the fall of potential is entirely due to the 
unscreened fluctuations of the disorder potential (Fig. 
||b). Note that while a purely electrostatic solution in the 
absence of disorder predicts a complete screening in the 
bulk for any filling factor, the situation is, in fact, con- 
siderably different due to the interplay of disorder and 
LL quantization. 

We conclude that at the left edge of the sample there is 
a geometric boundary line T which ought to be infinite in 
the y-direction, and which marks the transition from the 
last compressible strip to the incompressible bulk (Fig. 
2a). To the left of this line the density of the electron 
charge falls off, and there are N percolating compressible 
strips. We will refer to this region from the physical edge 
of the system to the boundary line T as to the 'edge chan- 
nel region' (there are two such regions - one at each edge 
as the argument above clearly holds for the right edge, 
as well). On the other side of this boundary an incom- 
pressible region with N completely filled LL's percolates. 
In this region the potential falls due to the unscreened 
fluctuations of the disorder potential (Fig. 2b). To this 
whole region we will refer to as the bulk. 

Let us first show that in the regime of the QHE no 
part of the non-equilibrium Hall current flows in the edge 




E(x) 




FIG. 2. Charge distribution and energies of Landau lev- 
els near the physical boundary a. Compressible and incom- 
pressible 'strips' near the edge. b. Energies of the first three 
occupied Landau levels close to the sample boundary. 



channel region as defined above (for a more detailed dis- 
cussion see Ref. 17). 

In equilibrium, although the net current is zero, orbital 
current flows along the trajectories extended from source 
to drain at one edge and in the opposite direction at the 
other edge. There are two contributions to—this equilib- 
rium current which are spatially separateco. First, the 
diamagnetic current, is roughly proportional to the gra- 
dient of the self-consistent confining potential, and it is 
non-zero only in the incompressible regions. The other 
one is associated with the concentration gradient, which 
is different from zero only in the compressible strips close 
to the device boundary. However, the total current car- 
ried by the partially occupied near-edge region of the n-th 
LL is always (n+l/2)ew c /27r, independent of the external 
conditional. Corresponding compressible strips at op- 
posite edges of the sample carry the same current but in 
opposite directions, therefore only incompressible regions 
can support non-equilibrium current. This can also be 
seen from classical considerations: non-equilibrium Hall 
current is carried by electrons drifting in crossed elec- 
tric and magnetic fields in the direction perpendicular to 
both of them. This drift velocity is proportional to the 
local electric field, which is identically zero in the com- 
pressible regions. We will, therefore, ignore the current 
in the compressible regions (even if extended trajectories 
are present there) since it is not changed by the injected 
non- equilibrium current. 

Consider now the case when current is injected into 
the sample. Under conditions of the QHE, the longi- 
tudinal conductivity is negligible, and there is no dis- 
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sipative current between the two edges. Therefore, for 
every value of the injected current smaller than the crit- 
ical current corresponding to the QHE breakdown, there 
exists a many-particle steady state that carries the cur- 
rent. This non-equilibrium steady state can be described 
by certain charge and current distributions in ways sim- 
ilar to an equilibrium state. The only difference is that 
in the steady state the net current in the transverse di- 
rection is different from zero. As a consequence of this, 
LL's are filled up to different energies at the left and the 
right edges, and a Hall voltage can be measured. As the 
non-equilibrium current is carried by the states in the 
incompressible regions, the total current flowing in the 
sample can be obtained by summing over the contribu- 
tions due to the electrons in the fully occupied extended 
states. 

Let us first compute the total amount of currept carried 
in the edge channel region. Following Trugmancil we cal- 
culate the Hall current carried by one LL through some 
horizontal cross section y — y of length AL = x-i — x\ 
along which the potential falls by AV — V2 — V\. The 
number of electrons per unit area n is equal to l/2irhl 2 = 
eB/hc. Define the expectation value of the local electron 
velocity operator at each point along the line x±xj?3: 



where £(r) is the electric field. The current I Xl x 2 t na t, 
passes through the line X1X2 is: 

I - rdxneiv a)= e2 n* a ,yb)-y(*i,yb) 

J Xl h e 

(3) 

Since local currents, circulating in the closed trajecto- 
ries, contribute neither to total diamagnetic current nor 
to non-equilibrium Hall current (we will study their im- 
portance in Section |v|), let us choose the points (xi,yo) 
and (x2,yo) on the extended trajectories at potentials 
V\ and V2 respectively. Then the strip between these 
two trajectories contributes the amount proportional to 
Ay = V2 — V\ to the total current as given by Eq. (||). 

The Hall conductivity is a xy = I/Vjj, where Vjj is the 
Hall voltage measured in the four-terminal resistance ex- 
periments. This voltage is determined by the amount of 
current carried by each LL just before the voltage con- 
tacts, as well as by the transmission and reflection coef- 
ficients of these contacts. It is possible to showcl that the 
Hall conductivity will be quantized only if, upon leav- 
ing the voltage contacts, all LL's are filled up to the 
same energy: e n> ; = /ii, e n>r — u 2 for all n. In this 
case, the voltmeter measures the difference between the 
chemical potentials at the two edges: Vh = ("2 — 
Then from Eq. (J3J) follows that a xy — Ne 2 /h. We must 
note at this point that filling of each LL starts when the 
energy of this LL drops to [i\ (/Z2) at the left (right) 
edge. This means that all partially occupied states at 



one edge have the same energy (/zi or u 2 ). Then, both 
in equilibrium and in the non-equilibrium steady state 
described above, the total drop of the potential AV n in 
the part of the edge channel region where the n-th LL is 
filled, equals to (N — n)huj c . Therefore, the potential falls 
by (n — \)fiu c before the n-th LL starts to fill, and by 
[N — l)fiLo c through the whole edge channel region. As 
the energy difference between the LL's in the compress- 
ible strips does not change when the current is injected, 
the total potential drop A V n is also not changed by the 
injection of the Hall current. It follows then from Eq. 
(||) that the current flowing in the edge channel region 
is always equal to the equilibrium diamagnetic current, 
even when the non-equilibrium current is injected into 
the sample. In other words, no non-equilibrium current 
flows in the edge channel region, and the entire injected 
current is carried by the states in the bulk of the system, 
to the right from the line F on Fig. ^a (see also Ref . [l7]) . 
The rest of the paper discusses the distribution of the 
extended trajectories in the bulk of the sample. 

It is important to realize that the complete equilibra- 
tion at the edge leads, therefore, to complete absence 
of non- equilibrium current in traditionally defined edge 
channels. On the other hand, if this equilibration is not 
perfect, some portion of the non- equilibrium current is 
pushed into the edge channels. However, in this case, the 
measured Hall conductance is not quantized in .units of 
e 2 /h. Experiments, involving non-ideal contactsE] (con- 
tacts designed to avoid complete equilibration) clearly 
demonstrated this absence of quantization and implied 
presence of currents in the edge channels. In fact, par- 
tial equilibration was still achieved in such experiments 
due to inter-channel scattering. It was argued, however, 
that the only soutce of perfect equilibration are the ideal 
voltage contactsO. We argue, additionally, that employ- 
ing the ideal contacts leads necessarily to the absence of 
non- equilibrium current in the edge channels. Therefore, 
the experimental conditions are crucial in determining 
the Hall current distribution. 



B. The extended states, a continuity theorem 

Here we determine the spatial location of the extended 
states under QHE conditions. Let us consider the behav- 
ior of the potential immediately to the right of the line T. 
As mentioned above, this potential initially falls due to 
the incompressibility in this region (see Fig. |^b) . Exactly 
how far the potential falls depends on the local config- 
uration of the random potential, and at some point the 
potential will necessarily start raising again. Therefore, 
there is a spatially extended region to the right of T where 
the potential falls. Since T itself is an equipotential line 
at energy er = /i — NTilj c , by continuity there is a finite 
range of energies ei ast < e < er in which equipotentials 
percolate in the y-direction. Here we denoted the lowest 
energy of these equipotentials by ei as t and the chemical 
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FIG. 3. Extended trajectories from the edge towards the 
bulk of the sample a. Open trajectories: F is the inner bound- 
ary of the edge region, F' is the last open trajectory in the 
bulk. b. Corresponding energy profiles along the marked 
cross section lines AA and BB. 

potential by [i. We now ask what determines this min- 
imum energy ei as t and the shape of the corresponding 
equipotential. The results are far from obvious. 

It was argued before that in the absence of the confin- 
ing effects of the physical boundary (we will call such an 
artificial construction a 'quasi- infinite system'), the only 
effect of the finite size in x-direction would be the exis- 
tence of a narrow band of equipotentials connecting left 
and right edges. Let us assume that an equipotential V 
with energy eo belongs to this band. It is useful to note 
that at this energy there may also be other trajectories 
that do not percolate. Going back to a more realistic sys- 
tem, with confining potential at the edges, one can see 
that this percolating trajectory T' is necessarily cut off 
by a rise of the potential towards the physical edge (Fig. 

. It is then easily seen that at each energy between eo 
and er there exists one and only one equipotential line 
open in y-direction, and that no trajectory at energies 
lower than eo can percolate vertically. 

Proof: Let us connect lines T and T' by means of an 
arbitrary line 7. Along this line the potential has to fall 
from er at T to eo at T'. The potential along the 7 line 
does not have to be monotonous. 

a. Let us initially assume that the potential is 
monotonous along 7. Now connect 7 with the left and 




FIG. 4. Illustration of the proof of continuity theorem. 

right edges by, say, horizontal straight lines AA\ and BB\ 
(Fig. ||). Choose any point P on 7 between A and B. We 
now ask whether an equipotential line threading through 
P percolates in the y-direction. In order for the trajec- 
tory passing through P to be a closed one, one has to 
find another point on the line A\B\ with the same po- 
tential as in point P. Such point cannot be found on 
7 because of the monotonous behavior of the potential 
along it. There could be such points on the piece BB\ 
but in order to reach them one has to cross line V which 
is at different potential. Hence, there is no way to close 
equipotential trajectory by going to the right of P. A 
similar argument shows that one cannot close trajectory 
to the left from P either. Therefore, the equipotential 
passing through any point on 7 is necessarily delocalized 
in the y-direction. 

b. Going back to the general case of a non-monotonous 
potential on 7, we argue that for any energy e such that 
eo < e < er there will be an odd number of points on 
7: namely, there will be one more point on the fall of 
potential than on the rise. Some of these points com- 
bine pairwise into closed orbits. However, the additional 
single point satisfies conditions similar to those stated 
in part a: the equipotential passing through it must be 
extended in the y-direction. 

From this discussion it is also clear that no trajectories 
with energies beyond the interval [eo, er] can percolate in 
^-direction. 

We immediately recognize that eo = ei ast and that the 
corresponding right-most vertically percolating equipo- 
tential is the line we denoted V on Fig. 3. Therefore, 
as a first approximation, the trajectory V has to perco- 
late in both vertical and horizontal directions. This is a 
critical trajectory with at least one saddle point O on it. 
An equilibrium picture of extended trajectories emerges 
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according to this discussion as shown in Fig. |3| 

Qualitatively, we conclude that at least the trajecto- 
ries with energies greater than, but sufficiently close to £q, 
travel throughout the system towards the 'saddle point', 
as does the critical trajectory itself. The latter comes 
near the boundaries only to get around closed trajec- 
tories with similar energies and with large sizes in the 
x-direction. Since such trajectories are very sparse, the 
critical trajectory spends most its path inside the system 
rather than near the boundary. In addition, there can 
also exist vertically percolating trajectories in the edge 
compressible regions to the left of line T (see Fig. ||) . 



C. General remarks 

Insofar, by means of purely geometrical arguments, we 
have proved that: 

i. Contrary to a widespread belief, it is not just the 
existence of the physical boundaries which is responsi- 
ble for the formation of a band of extended states (in- 
stead of a single extended state in the infinite random 
two-dimensional system). It is clear from the preceding 
discussion that the number, density and location of the 
trajectories extended in y-direction have little to do with 
the sample size in x-direction. This is in complete con- 
trast with the common statement that the bandwidth 
of the extended states is inversely proportional to the 
sample size. The most important factor is the interac- 
tion between electrons which leads to the screening of the 
confining and disorder potential at the physical bound- 
ary and, consequently, to the formation of an extended 
boundary line T. It is the difference between the energy 
of this equipotential line T (which is equal to the chemi- 
cal potential /i) and the energy of the state extended in a 
corresponding infinite system, that determines the total 
number and density of extended states in equilibrium. 

ii. In addition, we must conclude that even in equi- 
librium, extended states are present not only near the 
physical boundary but also deep in parts of the bulk, 
their particular, sample-dependent distribution being de- 
termined by a specific configuration of the disorder. 

In the next section we give a detailed description of 
how, under the conditions of complete equilibration at 
each edge, the Hall current is distributed in the sample. 



IV. DELOCALIZED TRAJECTORIES AND 
CURRENT DISTRIBUTION IN THE STEADY 
STATE 

Imagine going from the line T towards the inside of the 
system in the current carrying steady state. The differ- 
ence from the equilibrium state then is that the coun- 
terpart of r at the right edge is an equipotential at an 
energy lower than er by the Hall voltage Vh- Using ar- 
guments similar to the ones applied in equilibrium, it is 



easy to demonstrate that there is a critical trajectory V 
with at least one saddle point in the steady state. The 
energy of this critical trajectory will not be the same 
as in equilibrium: in our setting of Hall field directed 
from right to left, the following inequality should hold: 
(er-er')ixon-equil > (er-er')equil- Therefore, the number 
of extended trajectories differs from the one in equilib- 
rium. This should be true independently of whether the 
charge distribution is changed due to the presence of Hall 
electric field or not. 

As the Hall current is injected into the sample, it 
charges the edge metallic regions establishing the chem- 
ical potential difference between the two opposite edges. 
These charges produce the long-range electric field de- 
caying as the inveesepdjatance from each edge sufficiently 
far from this edgeuLESEl In the context of the problem 
of Hall current distribution the effecjUof charged edges 
was thoroughly studied in the pastoEl'E-X If these charges 
cannot move into the bulk they create a logarithmic 
potential sufficiently far from the edges. This solution 
assumes zero longitudinal conductivity. More realisti- 
cally, however, one has to estimate the charge relaxation 
times associated with the finite longitudinal conductiv- 
ity o~xx or with some other charge relaxation mechanisms 
(not related to the transport by the extended states) . 
ThoulessEj considered a finite a xx , and assumed as the 
condition for a steady state a constant value of dissipative 
current j x (x) — const. The problem reduces to the two- 
dimensional Poisson equation and is solved for the charge 
distribution inside the Hall bar. This new charge distri- 
bution produces a homogeneous electric field across the 
sample and can be established only if there are regions 
in the bulk which can absorb additional charges. The 
origin of such regions, their sizes and spatial distribution 
will be discussed at the end of the paper in the context 
of recent measurements, and for the estimates of typical 
charge relaxation times. A mo re .. thp c ough treatment can 
be found in other publicationscBEZIcj. 

Consider first the situation when the relaxation time 
is very long compared with the time of measurement. In 
this case the Hall electric field is negligible at the saddle 
point on the critical trajectory, falling off inversely with 
the distance from the edge. One can argue that this crit- 
ical trajectory changes its energy but the saddle point 
on it does not move. Consequently, all the newly formed 
percolating trajectories have to be located between the 
critical trajectory and the adjoining edge region, that is, 
in the region where delocalized states were present in 
equilibrium. Similarly to equilibrium delocalized trajec- 
tories, the newly formed ones are not distributed homoge- 
neously. While in most of the area of the sample their av- 
erage concentration follows the Hall electric field (by the 
arguments used to prove the equilibrium distribution), in 
the regions where the critical trajectory approaches the 
edge, the density of newly formed delocalized states has 
to be much higher. 

Consider a horizontal cross section AA\ where the crit- 
ical line r' approaches the line T at a distance A <C W 
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(Fig. 5). The fall of the potential AV on the distance A 
due to the external field is negligible compared to the Hall 
voltage. Yet, the application of a small Hall field cannot 
destroy many localized states inside the critical trajec- 
tory as in most parts of the system the localizing disorder 
field is much greater than the external one. These two 
facts can be reconciled if most of the newly formed de- 
localized trajectories with energies in the interval of the 
order of Vr appear in the narrow neighborhood of the 
original position of T', to the right of it. The trajectory 
r' itself shifts to the right by a distance 8 to let the verti- 
cally percolating trajectories penetrate into the interval 
AA\ with length 5 (Fig. ||). Such a solution is possible 
due to the fact that the slope of the potential in the inter- 
val AAi is determined by the field of disorder potential 
rather than by the much weaker external field. In the 
areas of the sample where the critical trajectory is suffi- 
ciently far from the edges so that the Hall field near this 
critical trajectory is negligible, the geometry of the crit- 
ical line does not affect the position of the field-induced 
delocalized states. The concentration of latter ones fol- 
lows the intensity of the Hall field and the concentration 
of equilibrium extended trajectories: new trajectories are 
formed as a response to the Hall electric field on those 
slopes of disorder potential which contained delocalized 
states in equilibrium. One concludes, therefore, that, 
except for the electrostatic Hall electric field, there are 
no specific features brought into the Hall current spatial 
distribution by the edges. Since equilibrium delocalized 
trajectories are generally not concentrated near the edge 
(previous Section), neither are the newly formed ones. 
The density of non-equilibrium extended trajectories is 
weighted towards the edges only to the extent the Hall 
electric field is. 

The other limiting case is the condition of complete 
(fast) relaxation. It is realized if the charges, including 
those on localized states, respond to the external elec- 
tric field created by the charged edges and redistribute 
themselves according to the steady state solution. The 
typical rate of this response must be greater than the rate 
at which the parameters are changed in the experiment. 
The only possible source of charge redistribution are the 
metallic regions of partially filled Landau levels. There- 
fore, one has to estimate the time it takes for a finite 
charge to be transfered from any such metallic region to 
the neighboring one. 

There is experimental evidence that, in many circum- 
stances, equilibration does occur. Direct study of relax- 
ation processesE3 shows that the state of an incompress- 
ible system achieved by raising the magnetic field from 
zero to several Tesla at high temperature (of the order of 
IK) differs significantly from the one reached by increas- 
ing the magnetic field in a pre-cooled (lOOmK) sample 
and then heating it-up back to IK. Although the relax- 
ation considered ir£3 requires spin-flip processes and is, 
therefore, significantly slower, one can think of the relax- 
ation as being facilitated by higher temperatures allowing 
for greater hopping or activated current. 
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FIG. 5. Critical trajectory F' approaching the boundary 
line F and the form of electrostatic potential in this region 
(see the text). 



Other group of experiments produce some indirect ev- 
idence of relaxation: the dependence of QH breakdown 
current on the sample width. Most of these measure- 
ments show a linear dependence which contradicts the 
assumption that the current flows primarily through the 
edges. This implies that that the electric field is close to a 
constant across the sample. This can be achieved only by 
delivering electrons to localized states in the compress- 
ible islands. At currents close to the critical current for 
breakdown of the QHE the Hall electric field is compara- 
ble with the field of disorderly and this factor itself can 
initiate effective relaxation processes via possible tunnel- 
ing into unoccupied states on the higher Landau levela. 
However, right before the breakdown, the bulk of the 
system is a well-defined incompressible system, in which 
relaxation processes are fast. We will present a consis- 
tent theoretical description and quantitative conclusions 
in the next Section. 

Consider now where the delocalized trajectories are lo- 
cated in such an equilibrated system in a weak Hall elec- 
tric field. The picture of extended trajectories is changed 
relative to the situation of an insulating (rigid) bulk only 
by the fact that the electric field is homogeneous on av- 
erage rather than weighted towards the edges. This dif- 
ference, however, leads to the distribution of extended 
states being on average a copy of such a distribution in 
equilibrium. Therefore, the current injected into the sys- 
tem in such case flows in the region between the lines T 
and r' as they were defined in equilibrium and has no 
preference of the edge regions. In the other words, the 
current is distributed throughout the system. 
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V. IMPLICATIONS OF CURRENT 
DISTRIBUTION 

In either of the two limits of short or long relaxation 
times, for sufficiently small Hall currents, one can label 
the corresponding regime as 'linear response'. In such a 
'linear response' regime the critical trajectory V plays 
the role attributed in the edge-channels formalism to the 
innermost boundary of the edge channels. However, un- 
like the edge channels, V itself is not confined to be near 
the physical boundaries. In this sense the description 
proposed here is complementary to the edge state for- 
malism. It also resolves _spme controversy around the 
validity of Stfeda formulao. 

In the derivation of the Stfeda formula it is assumed 
that the derivative (dn/dB)^ is quantized due to the ex- 
istence of the gap in the density of states between Lan- 
dau levels. In this formula n stands for the density of 
states, B for magnetic field, and /i for the fixed chemical 
potential. This assumption is crucial for the use of the 
Maxwell relation (dn/dB)^ = (l/A)(dM/dfj,) B , where 
M is the magnetization and A is the total area of the 
system. Under certain assumptions (which we mention 
below) the derivative in the right-hand side is equal to 
conductivity thus providing the proof for the conductiv- 
ity quantization. It is clear, however, that there is just 
a mobility gap between LL's but there is a finite den- 
sity of localized states in this mobility gap. These local- 
ized states have to contribute to the derivative (dn/dB)^ 
in one way or another. This immediately questions the 
validity of the assumptions used to prove the equality 
o~ X y = (l/A)(dAl/d[i)B , namely that the change in mag- 
netization SM is entirely due to additional edge current 
associated with the change in chemical potential 8fi. 

For interacting electrons, however, 8M is constructed 
from both drift along extended equipotentials and circu- 
lating currents on closed loops. The latter ones do not 
contribute to the Hall conductivity but turn out to be 
essential in evaluating-both right- and left-hand sides of 
the Maxwell relatiorO. Our argument makes use of the 
existence of states in the mobility gap, on the one hand 
by implying the formation of 'new' extended trajecto- 
ries within this gap and on the other, by explaining how 
these new extended trajectories become responsible for 
the total change in magnetization. 

At higher Hall voltages one expects the original topol- 
ogy of states to be gradually destroyed by the external 
field. The profile of the random potential has to change to 
absorb more delocalized trajectories than it would be able 
to absorb in any of the two scenarios described above. In 
fact, only at rather strong Hall currents does the treat- 
ment presented in this Section need to be reconsidered. It 
happens only when the relaxation processes and the dis- 
sipative transport change significantly compared to the 
linear response regime. This, in turn, occurs when the 
local Hall electric field, at least in some areas of the sam- 
ple, becomes comparable with the electric field of the 



disorder potential. This field turns out to be of the order 
of 0.1Tllu c /Ib. The corresponding to Hall current den- 
sities are lOA/m which exceeds experimental values in 
most cases. At such strong external fields the system 
becomes more homogeneous as the potential fluctuations 
get wiped out. The newly formed percolating trajecto- 
ries spread away from the critical line T' and distribute 
throughout the sample. 



VI. INTERPRETATION OF RECENT 
EXPERIMENTS 

A theoretical model which captures most of the quali- 
tative and quantitative features of charge relaxation pro- 
cesses in the QH regime has been constructed in Refs. 

,0. Here we will sketch the main blocks of the model 
and present the results leading to.-a.-cpnsistent explana- 
tion of the two recent experimentsoli3. Here we concen- 
trate on the experiments of Ref. |l^. 

The random resistor-capacitor network model 
(RRCNM) is based on the theory of charge transfer be- 
tween isolated metallic islands in the bulk via hopping. 
Due to the small size of these islands one has to take 
into account the Coulomb blockade effects: adding or 
removing a single electron from the island costs a non- 
negligible charging energy. Therefore, each island serves 
as a capacitor in the model. Typically the charging en- 
ergies are of the order of 1 meV. Moving electron from 
a charged island to an adjacent neutral one leads to ei- 
ther emission or absorption of a phonon with the energy 
equal to the difference of the charging energies of the 
two islands. This difference is typically 0.3 — 0.5 meV, 
or 3 — 5 K. Such phonon-assisted hopping is somewhat 
suppressed, in a case of phonon absorption. However, in 
the range of temperatures between 300 mK and 3 K in 
which most of the experiments on current distribution 
and charge relaxation are performed, this suppression is 
not dramatic. 

The other factor limiting the rate of charge transfer 
through the bulk is the distance between the communi- 
cating islands. The extent of the wave function of an 
electron at the boundary of a compressible island is of 
the order of a few magnetic lengths. If the distance be- 
tween the islands significantly exceeds this width, the 
overlap between the states participating in an individual 
act of hopping is exponentially small leading to a van- 
ishing hopping probability. This limitation defines a few 
(if any) paths which are available for a charge transfer. 
Note that this process not only leads to a finite longitu- 
dinal current but also to the formation of a new steady 
state with a charge distribution different from the one in 
equilibrium. The links between the islands enter into the 
model as resistive elements. 

The average distance between adjacent islands varies 
substantially with the filling factor. In the range of tem- 
peratures mentioned above, the QH plateau is quite nar- 
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FIG. 6. The basic block of the random resistor-capacitor 
network. 



row: it spans just 5% to 20% of the entire distance be- 
tween the Landau levels on the scale of filling factors. 
The greater the number of filled LL's, the narrower the 
plateau and the smaller the typical distance between the 
compressible islands. It turns out that relaxation can 
be very slow for the filling factor v — 2 (spin-unresolved 
LL's) at T ~ 300 mK. However, for v — 4, even at these 
low temperatures the RRCNM predicts relaxation times 
of the order of 1 ms which is shorter than the typical 
experimental time scales. 

The basic block of the simplest RRCN is shown on 
Fig. [| Each node corresponds to a compressible island 
which is connected resistively to a neighboring one. The 
capacitance of each island to the ground is accounted for 
by the capacitors Ci while mutual capacitances between 
the islands i and j are denoted as Cy . 

The RRCNM is solved numerically for each particu- 
lar configuration of disorder and at different filling fac- 
tors. The solution provides the geometric paths of elec- 
tron transfer, the charges on compressible islands, AC 
and DC responses of the system, and charge relaxation 
time t. Given the conditions of an experiment, one can 
obtain these characteristics to predict or explain the ob- 
served results. j— - 

In recent measurementsEj, the atomic force microscope 
technique has been used to study the distribution of the 
Hall current in the QH sample. Their results which we 
interpret below, are shown on Fig. |t]. 

The scheme of this e xperim ent is different from the 

This is not a transport- 



one discussed in section [II A 



type measurement which means that there are no volt- 
age contacts attached to the sample. The voltage profile 
is measured capacitively, without a direct contact with 
the 2DEG. Therefore, the arguments of the section [II A 
proving the absence of non-equilibrium Hall current in 
the edge channels, will not apply to this experiment in 
general. More specifically, there is no ideal equilibra- 
tion between the edge channels provided by the voltage 
contacts in transport measurements. The only source of 
inter-edge-abannel equilibration is the mutual scattering 
of electrons^. It is less effective than the equilibration 
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FIG. 7. Voltage profiles across a 12 /im-wide Hall bar on a 
low-mobility sample near the v = 6 QH plateau (traces offset 
for clarity). The B- fields at which the profiles were measured, 
and the v — 6 plateau measured by transport are shown in the 
inset, (a) Well below the plateau, the Hall voltage profile is 
linear. (b,c) As the plateau is approached, the bulk retains a 
linear profile but decouples from the edges, (d-g) Well inside 
the plateau, the profile in the bulk becomes complicated and 
changes rapidly with field. (h,i) A linear profile returns upon 
leaving the plateau, [after Ref. [j^ 



in the voltage contacts but for sufficiently narrow edge 
channels it still leads to the establishment of equal chemi- 
cal potentials in the compressible edge regions of different 
LL's. 

On Fig. ^ the results of Hall resistance R xy measure- 
ments are plotted. The curve ^(a) corresponds to a 
filling factor v ~ 6.8. The bulk is clearly a metallic 
system. There are three (spin-degenerate) compressible 
edge channels (0 < v < 2, 2 < v < 4, and 4 < v < 6) 
and three incompressible edge channels (v— 2,4,6). As 
shown in Ref. [3l|, only the inner-most edge regions can 
be wide - up to 1§Ib- Since the resolution of the exper- 
iment is about 2QIb : the edge regions, whether they are 
equilibrated or not, cannot be resolved. This leads to a 
linear profile of the Hall voltage. 

As the QH plateau v = 6 is approached with the in- 
creasing magnetic field, the bulk filling factor decreases. 
On the curve f?](b) the bulk filling factor is vq ~ 6.2. 
The filling profile is sketched on Fig. [| One can use 
the electrostatic solution for the 2DEG density at the 
edge, Eq. (m), to estimate how fast the top-most 7-th 
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FIG. 8. Interpretation of the experiment in Fig. @ (Ref. 
|l^). The structure of the 2DEG near the edge. Top: al- 
ternating compressible (shaded) and incompressible (blank) 
regions. Bottom: filling factor profile. The dashed line on 
the top picture and point 2 on the bottom correspond to the 
distance where the top-most-S-th LL starts filling according 
to the electrostatic solutionEil. In reality, the preceding in- 
compressible edge channel with v = 6 extends deeper into the 
bulk (point 3) due to the presence of disorder (see the text). 



LL gets filled to its bulk value. Assuming for the de- 
pletion width I the typical value I ~ 101b, we find that, 
for vq = 6.2 the distance between the line with v = 6.0 
(point 2 on Fig. ||) and v — 6.1 (point 3 on the same pic- 
ture) is approximately 3001b compared to just 101b for 
Vq = 6.8. Both experimental evidence of the filial widths 
of QH plateaus and our numerical simulationscJ clearly 
demonstrate that the system remains incompressible in 
a certain range of filling factors Av around the integer. 
Given the parameters used in these experiments (mobil- 
ity 80000 cm 2 /Vs, density 2.8 xlO 11 cm~ 2 and temper- 
atures between 0.7-1.0 K), according to our numerical 
simulations (details of which can be found in Ref. 29), 
Av > 0.1. Therefore, on a length of the order of 3001b 
to the right from point 2 on Fig. || the system is still 
incompressible. Such a wide incompressible strip makes 
mutual equilibration between the inner-most compress- 
ible edge channel (4 < v < 6) and the outer-most states 
of the compressible bulk practically impossible. To find 
out where the Hall voltage accumulates one has to study 
local longitudinal conductivity a xx . Since the bulk is 
metallic, its a xx is quite large (typical metallic conduc- 
tivity is of the order of 0.1e 2 /h). On the other hand, the 
longitudinal conductance of a wide incompressible strip 
is exponentially small. The Hall current flows in the re- 
gion with smaller dissipative conductivity. This explains 
the observation on Fig. ^(b) , (c) of large potential drops 
near the edges. 



As the bulk filling factor approaches vq = 6 even closer, 
the bulk becomes less and less conductive and the incom- 
pressible edge region gets wider (eventually connecting 
the two edges and transforming the bulk into the incom- 
pressible state at Vo ss 6.1). This leads to a smooth 
transition from an almost linear potential profile across 
the entire sample Fig. |?](d),(e) to a regime where the 
potential is flat at the edges and falls only in the bulk 
(Fig. |?](f),(g)). The profile of the potential in the last 
regime confirms two conclusions we have made earlier in 
this paper. The linear fall of the potential in the bulk 
means there is a complete charge relaxation. In fact, the 
numerical calculations based on the RRCNM predict re- 
laxation times of the order of 10 _5 s for the experimental 
conditions. This is by far the shortest time present in 
the problem. Second, oscillations of the curves around 
the linear profile are caused by the inhomogeneities in 
the bulk: a linear average electrostatic potential is pro- 
duced by the charged metallic islands which vary in size 
and in separation from each other. 



VII. CONCLUSIONS 



We have proved, by means of purely geometrical ar- 
guments, that a band of extended states, extended from 
source to drain, results from the interplay between elec- 
tron interactions and random and confining potentials. 
These delocalize d stat es extend through the bulk of the 
system (see Sec. IIIC). 

Based on these conclusions we built a consistent de- 
scription of the delocalized states in the regime when 
classical percolation theory is applicable. We were able 
to demonstrate that, under the most general conditions, 
the Hall current-carrying states are distributed in the 
bulk of the Hall sample rather than confined to the edge 
region. We showed how particular experimental condi- 
tions influence the results of the measurements, and yet 
confirm the general picture of distributed Hall currents. 
We believe that the controversy regarding the Hall cur- 
rent distribution, as well as the nature and location of 
extended states in the quantum Hall system is now re- 
solved. 



The low temperature regime, however, still needs to 
be addressed since the equivalence between equipotcn- 
tials and electron trajectories is lost when quantum ef- 
fects become essential (for an extended discussion on the 
transition between the two regimes see Refs. |29|j42] ). The 
authors do not know of any experiment on the distribu- 
tion of the Hall current and electric fields which explores 
this low temperature regime. Theoretically a more com- 
plete understanding of the microscopic structure of dis- 
ordered quantum Hall liquid is necessary to solve this 
problem. 
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